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We present new full-sky temperature and polarization maps in five frequency 
bands from 23 to 94 GHz, based on data from the first five years of the WMAP sky 
survey. The new maps are consistent with previous maps and are more sensitive. 
The five-year maps incorporate several improvements in data processing made 
possible by the additional years of data and by a more complete analysis of 
the instrument calibration and in-flight beam response. We present several new 
tests for systematic errors in the polarization data and conclude that W band 
polarization data is not yet suitable for cosmological studies, but we suggest 
directions for further study. We do find that Ka band data is suitable for use; 
in conjunction with the additional years of data, the addition of Ka band to the 
previously used Q and V band channels significantly reduces the uncertainty in 
the optical depth parameter, r. Further scientific results from the five year data 
analysis are presented in six companion papers and are summarized in §7] of this 
paper. 

With the 5 year WMAP data, we detect no convincing deviations from 
the minimal 6-parameter ACDM model: a flat universe dominated by a cos- 
mological constant, with adiabatic and nearly scale- invariant Gaussian fluctua- 
tions. Using WMAP data combined with measurements of Type la supernovae 
(SN) and Baryon Acoustic Oscillations (BAO) in the galaxy distribution, we 
find (68% CL uncertainties): tt b h 2 = 0.02267^;S, tt c h 2 = 0.1131 ± 0.0034, 
tt A = 0.726 ± 0.015, n s = 0.960 ± 0.013, r = 0.084 ± 0.016, and = (2.445 ± 
0.096) x 10" 9 at k = 0.002 Mpc" 1 . From these we derive: a 8 = 0.812 ± 0.026, 
H = 70.5 ± 1.3 km s" 1 Mpc" 1 , Q b = 0.0456 ± 0.0015, Q c = 0.228 ± 0.013, 
tt m h 2 = 0.1358tHo36; -Sreion = 10.9 ± 1.4, and t = 13.72 ± 0.12 Gyr. The new 
limit on the tensor-to-scalar ratio is r < 0.22 (95% CL), while the evidence for a 
running spectral index is insignificant, dn s /dlnk = — 0.028±0.020 (68% CL). We 
obtain tight, simultaneous limits on the (constant) dark energy equation of state 
and the spatial curvature of the universe: —0.14 < 1 + w < 0.12 (95% CL) and 
-0.0179 < fife < 0.0081 (95% CL). The number of relativistic degrees of free- 
dom, expressed in units of the effective number of neutrino species, is found to 
be A/gff = 4.4 ± 1.5 (68% CL), consistent with the standard value of 3.04. Models 
with N e g = are disfavored at >99.5% confidence. Finally, new limits on phys- 
ically motivated primordial non-Gaussianity parameters are —9 < /]^£ al < 111 
(95% CL) and -151 < f^f < 253 (95% CL) for the local and equilateral models, 
respectively. 

Subject headings: cosmic microwave background, cosmology: observations, early 
universe, dark matter, space vehicles, space vehicles: instruments, instrumenta- 
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tion: detectors, telescopes 

1. INTRODUCTION 

The Wilkinson Microwave Anisotropy Probe (WMAP) is a Medium-Class Explorer 
(MIDEX) satellite aimed at elucidating cosmology through full-sky observations of the cosmic 
microwave background (CMB). The WMAP full-sky maps of the temperature and polariza- 
tion anisotropy in five frequency bands provide our most accurate view to date of conditions 
in the early universe. The multi-frequency data facilitate the separation of the CMB signal 
from foreground emission arising both from our Galaxy and from extragalactic sources. The 
CMB angular power spectrum derived from these maps exhibits a highly coherent acoustic 
peak structure which makes it possible to extract a wealth of information about the compo- 
sition and history of the universe, as well as the processes that seeded the fluctuations. 
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model of cosmology: a flat universe dominated by dark energy, supplemented by dark matter 
and atoms with density fluctuations seeded by a Gaussian, adiabatic, nearly scale invariant 
process. The basic properties of this universe are determined by five numbers: the density of 
matter, the density of atoms, the age of the universe (or equivalently, the Hubble constant 
today), the amplitude of the initial fluctuations, and their scale dependence. 



By accurately measuring the first few peaks in the angular power spectrum and the 
large-scale polarization anisotropy, WMAP data have enabled the following inferences: 

• A precise (3%) determination of the density of atoms in the universe. The agreement 
between the atomic density derived from WMAP and the density inferred from the 
deuterium abundance is an important test of the standard big bang model. 

• A precise (3%) determination of the dark matter density. (With five years of data 
and a better determination of our beam response, this measurement has improved 
significantly.) Previous CMB measurements have shown that the dark matter must 
be non-baryonic and interact only weakly with atoms and radiation. The WMAP 
measurement of the density puts important constraints on supersymmetric dark matter 
models and on the properties of other dark matter candidates. 
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A definitive determination of the acoustic scale at redshift z = 1090. Similarly, the re- 
cent measurement of b aryon acoustic oscillations (BAO) in the galaxy power spectrum 
( lEisenstein et al.l 120051 ) has determined the acoustic scale at redshift z ~ 0.35. When 
combined, these standard rulers accurately measure the geometry of the universe and 
the properties of the dark energy. These data require a nearly flat universe dominated 
by dark energy consistent with a cosmological constant. 

A precise determination of the Hubble Constant, in conjunction with BAO observa- 
tions. Even when allowing curvature (Qq ^ 1) and a free dark energy equation of 
state (w 7^ —1), the acoustic data determine the Hubble constant to within 3%. The 
measured val ue is in excellent agre ement with independent results from the Hubble 
Key Project (IFreedman et al.l 120011 ). providing yet another important consistency test 
for the standard model. 

Significant constraint of the basic properties of the primordial fluctuations. The anti- 
correlation seen in the temperature/polarization (TE) correlation spectrum on 4° scales 
implies that the fluctuations are primarily adiabatic and rule out defect m odels and 
isocurvature models as the primary source of fluctuations (jPeiris et al.l 120031 ) . 



Further, the WMAP measurement of the primordial power spectrum of matter fluctua- 
tions constrains the physics of inflation, our best model for the origin of these fluctuations. 
Specifically, the 5 year data provide the best measurement to date of the scalar spectrum's 
amplitude and slope, and place the most stringent limits to date on the amplitude of tensor 
fluctuations. However, it should be noted that these constraints assume a smooth function of 
scale, k. Certain models with localized structure in P(k), and hence additional parameters, 
are no t ruled out, but neither are they required by the data; see e.g. IShafieloo fc Souradeep 
(l2007h : iHunt fc Sarkari J2007h . 



The statistical properties of the CMB fluctuations measured by WMAP are close to 
Gaussian; however, t here are several h ints of possible d e viations from Gaussianity, e.g . 
Eriksen et aD J2007ah : Icopi et al.l ; lLand fc Magueiiol fl2~0~^ : lYadav fc WandeltN2008h . 

Significant deviations would be a very important signature of new physics in the early uni- 
verse. 

Large-angular-scale polarization measurements currently provide our best window into 
the universe at z ~ 10. The WMAP data imply that the universe was reionized long 
before the epoch of the oldest known quasars. By accurately constraining the optical depth 
of the universe, WMAP not only constrains the age of the first stars but also determines 
the amplitude of primordial fluctuations to better than 3%. This result is important for 
constraining the growth rate of structure. 
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This paper summarizes results compiled from 5 years of WMAP data that are fully 
presented in a suite of 7 papers (including this one). The new results improve upon pre- 
vious results in many ways: additional data reduces the random noise, which is especially 
important for studying the temperature signal on small angular scales and the polarization 
signal on large angular scales; five independent years of data enable comparisons and null 
tests that were not previously possible; the instrument calibration and beam response have 
been much better characterized, due in part to improved analyses and to additional years of 
data; and, other cosmological data have become available. 

In addition to summarizing the other papers, this paper reports on changes in the 
WMAP data processing pipeline, presents the 5 year temperature and polarization maps, 
and gives new res ults on in s trume nt calibration and on potential systematic errors in the 
polarization data. iHill et al.l (120081 ) discuss the program to derive an improved physical optics 
model of t he WMAP t elesco pe, and use the results to better determine the WMAP beam 



response. iGold et al.l (120081 ) present a new analysis of diffus e foreground emissio n in the 
WMAP data and update previous analyses using 5 year data. IWright et al.l (120081 ) analyze 
extragalact ic point sources an d provide an updated source catalog, with new results on source 
variability. iNolta et al.l (120081 ) derive th e angular powe r spect ra from the maps, including the 
TT, TE, TB, EE, EB, and BB spectra. iDunkley et al.l ( 120081 ) produce an updated likelihood 
function and present cosmological parameter results based on 5 year WMAP data. They 
also develop an independent analysis of polarized foregrounds and use thos e results to test 



the re liability of the optical depth inference to foreground removal errors. iKomatsu et al. 
( 120081 ) infer cosmological parameters by combining 5 year WMAP data with a host of other 
cosmological data and discuss the implications of the results. Concurrent with the submission 
of these papers, all 5 year WMAP data are made available to the research community via 
NASA's Legacy Archive for Microwave Background Data Analysis (L AMBDA). The dat a 
products are described in detail in the WMAP Explanatory Supplement (ILimon et al.ll2008l ). 
also available on LAMBDA. 

The WMAP instrument is composed of 10 diff erencing; assemblies (DAs) spanning 5 
frequencies from 23 to 94 GHz faennett et alil2003h : 1 DA each at 23 GHz (Kl) and 33 
GHz (Kal), 2 each at 41 GHz (Q1,Q2) and 61 GHz (V1,V2), and 4 at 94 GHz (W1-W4). 
Each DA is formed from two differential radiometers which are sensitive to orthogonal linear 
polarization modes; the radiometers are designated 1 or 2 (e.g., VI 1 or W12) depending on 
polarization mode. 

In this paper we follow the notation convention that flux density is S ~ v a and antenna 
temperature is T ~ z/^ 3 , where the spectral indices are related by (3 = a — 2. In general, the 
CMB is expressed in terms of thermodynamic temperature, while Galactic and extragalactic 
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foregrounds are expressed in antenna temperature. Thermodynamic temperature differences 
are given by AT = ATj\[(e x — l) 2 /x 2 e x ], where x = hu/kTo, h is the Planck constant, v 
i s the frequency, k is the Boltzmann constant, and T = 2.725 K is the CMB temperature 
(IMather et al.lll999l ). A WMAP band-by-band tabulation of the conversion factors between 
thermodynamic and antenna temperature is given in Table [TJ 



2. CHANGES IN THE 5 YEAR DATA ANALYSIS 

The 1 year and 3 year data analyses were described in detail in previous papers. In large 
part, the 5 year analysis employs the same methods, so we do not repeat a detailed processing 
description here. However, we have made several improvements that are summarized here 
and described in more detail later in this paper and in a series of companion papers, as 
noted. We list the changes in the order they appear in the processing pipeline: 

• There is a ~ 1' temperature-dependent pointing offset between the star tracker coordi- 
nate system (which defines spacecraft coordinates) and the instrument boresights. In 
the 3 year analysis we introduced a correction to account for the elevation change of 
the instrument boresights in spacecraft coordinates. With additional years of data, we 
have been able to refine our thermal model of the pointing offset, so we now include 
a small (<1') correction to account for the azimuth change of the instrument bore- 
sights. Details of t he new correction are given in the 5 year Explanatory Supplement 



flLimon et all 120081 ). 



We have critically re-examined the relative and absolute intensity calibration proce- 
dures, paying special attention to the absolute gain recovery obtainable from the mod- 
ulation of the CMB dipole due to WMAP's motion. We describe the revised procedure 
in §H and note that the sky map calibration uncertainty has decreased from 0.5% to 
0.2%. 

The WMAP beam response has now been measured in 10 independent "seasons" of 
Jupiter observations. In the highest resolution W band channels, these measurements 
now probe the beam response ~44 dB down from the beam peak. However, there is still 
non-negligible beam solid angle below this level (~0.5%) that needs to be measured to 
enable accurate cosmological inference. In the 3 year analysis we produced a physical 
optics model of the A-side beam response starting with a pre-flight model and fitting 
in-flight mirror distortions to the flight Jupiter data. In the 5 year analysis we have 
extended the model to the B-side optics and, for both sides, we have extended the 
fit to include distortion modes a factor of 2 smaller in linear scale (4 times as many 
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modes). The model is used to augm ent the flig h t bea m maps below a given threshold. 
The details of this work are given in iHill et al.l (120081 ). 



The far sidelobe response of the beam was determined from a co mbination of groun d 
measurements and in-flight lunar data taken early in the mission (IBarnes et al.ll2003l ). 
For the current analysis, we have replaced a small fraction of the far sidelobe data with 
the physical optics model described above. We have also m ade the following changes 
in our handling of the far sidelobe pickup (jHill et al.l 120081 ) : 1) We have enlarged the 
"transition radius" that defines the boundary between the main beam and the far 
sidelobe response. This places a larger fraction of the total beam solid angle in the 
main beam where uncertainties are easier to quantify and propagate into the angular 
power spectra. 2) We have moved the far sidelobe deconvolution into the combined 
calibration and sky map solver (Q. This produces a self-consistent estimate of the 
intensity calibration and the deconvolved sky map. The calibrated time-ordered data 
archive has had an estimate of the far sidelobe response subtracted from each datum 
(as it had in the 3 year processing). 

We have updated the optimal filters use d in the final step o f map-making. The func- 
tional form of the filter is unchanged (IJarosik et al.l 120071 ). but the fits have been 
updated to cover years 4 and 5 of the flight data. 

Each WMAP differencing assembly consists of two radiometers that are sensitive to 
orthogonal linear polarization states. The sum and difference of the two radiome- 
ter channels split the signal into intensity and polarization components, respectively. 
However, the noise levels in the two radiometers are not equal, in general, so more 
optimal sky map estimation is possible in theory, at the cost of mixing intensity and 
polarization components in the process. For the current analysis, we investigated one 
such weighted algorithm and found that the polarization maps were subject to unac- 
ceptable contamination by the intensity signal in cases where the beam response was 
non-circular and the gradient of the intensity signal was large, e.g., in K band. As a 
result, we reverted to the unweighted (and unbiased) estimator used in previous work. 

We have improved the sky masks used to reject foreground contamination. In previous 
work, we defined masks based on contours of the K band data. In the 5 year analysis 
we produce masks based jointly on K band and Q band contours. For a given sky cut 
fraction, the new masks exclude flat spect rum (e.g. free-free ) emission more effectively. 
The new masks are described in detail in iGold et al.l (120081 ) and are provided with the 
5 year data release. In addition, we have modified the "processing" mask used to 
exclude very bright sources during sky map estima tion. The new ma sk is defined in 
terms of low-resolution (r4) HEALPix sky pixels (IGorski et al.l 120051 ) to facilitate a 



- 8- 



cleaner definition of the pixel-pixel inverse covariance matrices, N^ 1 . One side effect 
of this change is to introduce a few r4-sized holes around the brightest radio sources 
in the analysis mask, which incorporates the processing mask as a subset. 



We have amended our foreground analysis in the following ways: 1) iGold et al.l (120081 ) 
perform a pixel-by-pixel analysis of the joint temperature and polarization data to 
study the breakdown of the Galactic emission into physical components. 2) We have 
up dated some aspec ts of t he Maximum Entropy (MEM) based analysis, as described 



m l Gold et al.l (120081 ). 3) iDunkley et al.l (120081 ) develop a new analysis of polarized 
foreground emission using a Gibbs sampling approach th at yields a cle a ned C MB po- 
larization map and an associated covariance matrix. 4) IWright et al.l (120081 ) update 
the WMAP point source catalog and present some results on variable sources in the 5 
year data. However, the basic cosmological results are still based on maps that were 
cleaned with the same template-based procedure that was used in the 3 year analysis. 

We have improved the final temperature pow er spectrum, CT T , b y using a Gibbs- 



based maximum likeli hood estimate for I < 32 (IDunkley et al.ll2008l ) and a pseudo-C; 
estimate for higher I (INolta et al.l 120081 ) . As with the 3 year analysis, the pseudo-C; 
estimate uses only V- and W-band data. With 5 individual years of data and six V- 
and W-band differencing assemblies, we can now form individual cross-power spectra 
from 15 DA pairs within each of 5 years and from 36 DA pairs across 10 year pairs, 
for a total of 435 independent cross-power spectra. 

In the 3 year analysis we developed a pseudo-C/ method for evaluating polarization 
power spectra in the presence of correlated noise. In the present analysis we addition- 
ally estimate the TE, TB, EE, EB, & BB spectra and th eir errors using an extension 
of the maximum likelihood method in iPage et al.l (120071 ). However, as in the 3 year 
analysis, the likelihood of a given model is still evaluated directly from the polarization 
maps using a pixel-based likelihood. 

We have improved the form of the likelihood f unction used to i nfer c osmological param- 

We use an exact 



eters from the Monte Carlo Markov Chains (IDunkley et al 



2008 



maximum likelihood form for the I < 32 TT data (jEriksen et al.l l2007ct ) . We have 
investigated theoretically optimal methods for incorporating window function uncer- 
tainties into the likelihood, but in tests with simulated data we have found them to 



be bi ased. In the end, we adopt the form used in the 3 year analysis (IHins 



law et al. 



2007), but we incorporate the smaller 5 year window function uncertainties (jHill et al. 
20081 ) as inputs. We now routinely account for gravitational lensing when assessing 
parameters, and we have added an option to use low-/ TB and EB data for testing 
non-standard cosmological models. 
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For testing nongaussiani ty, we employ an improved estimator for f^L (ICreminelli et al. 



2006 



fl2008h . 



Yadav et al.l 120071 ). The results of this analysis are described in 



Komatsu et al. 



OBSERVATIONS AND MAPS 



The 5 year WMAP data encompass the period from 00:00:00 UT, 10 August 2001 (day 
number 222) to 00:00:00 UT, 9 August 2006 (day number 222). The observing efficiency 
during this time is roughly 99%; Table [2] lists the fraction of data that was lost or rejected 
as unusable. The Table also gives the fraction of data that is flagged due to potential 
contamination by thermal emission from Mars, Jupiter, Saturn, Uranus, and Neptune. These 



data are not used in map-making, but are useful for in-flight beam mapping (jHill et al.ll2008 
Limon et al1l2008h . 



After performing an end-to-end analysis of the instrument calibration, si ngle-year sky 



maps are created from the time-ordered data using the procedure described by lJarosik et al. 
(120071 ). Figure [1] shows the 5 year temperature maps at each of the five WMAP observing 



frequencies: 23, 33, 41, 61, and 94 GHz. The number of independent observations per pixel, 
N ohs , is qualitatively the same as Figure 2 of iHinshaw et al.l (120071 ) and is not reproduced 
here. The noise per pixel, p, is given by a(p) = cr N ohs ' (p), where a is the noise per 
observation, given in Table HI To a very good approximation, the noise per pixel in the 5 
year maps is a factor of Vo times lower than in the single- year maps. Figures [2] and [3] show 
the 5 year polarization maps in the form of the Stokes parameters Q and U, respectively. 
Maps of the relative polar iz ation sensitivity, the Q and U analogs of A b s , are shown in 
Figure 13 of IJarosik et al.l (120071 ) and are not updated here. A description of the low- 
resoluti on pixel-pixel invers e covariance matrices used in the polarization analysis is also 



given in IJarosik et al.l (120071 ). and is not repeated here. The polarization maps are dominated 
by foreground emission, primarily synchrotron emission from the Milky Way. Figure [4] shows 
the polarization maps in a form in which the color scale represents polarized intensity, P = 
\J Q 2 + U 2 , and the line segments indicate polarization direction for pixels with a signal-to- 
noise ratio greater than 1. As with the temperature maps, the noise per pixel in the 5 year 
polarization maps is Vo times lower than in the single-year maps. 

Figure [5] shows the difference between the 5 year temperature maps and the correspond- 
ing 3 year maps. All maps have been smoothed to 2° resolution to minimize the noise 
difference between them (due to the additional years of data). The left column shows the 
difference without any further processing, save for the subtraction of a relative offset between 
the maps. Table [3] gives the value of the relative offset in each band. Recall that WMAP is 
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insensitive to absolute temperature, so we adopt a convention that sets the zero level in each 
map based on a model of the foreground emission at the galactic poles. While we have not 
changed conventions, our 3 year estimate was erroneous due to the use of a preliminary CMB 
signal map at the time the estimate was made. This error d id not affect any cosm ological 
results, but it probably explains the offset differences noted by lEriksen et al.l (j2007bl ) in their 
recent analysis of the 3 year data. 

The dominant structure in the left column of Figure consists of a residual dipole and 
galactic plane emission. This reflects the updated 5 year calibration which has produced 
changes in the gain of order 0.3% compared to the 3 year gain estimate (see §H]for a more 
detailed discussion of the calibration). Table [3] gives the dipole amplitude difference in 
each band, along with the much smaller quadrupole and octupole power difference. (For 
comparison, we estimate the CMB power at I — 2, 3 to be 1(1 + l)Q/27r = 211, 1041 /iK 2 , 
respectively.) The right column of Figure [5] shows the corresponding sky map differences 
after the 3 year map has been rescaled by a single factor (in each band) to account for 
the mean gain change between the 3 and 5 year calibration determinations. The residual 
galactic plane structure in these maps is less than 0.2% of the nominal signal in Q band, 
and less than 0.1% in all the other bands. The large scale structure in the band-averaged 
temperature maps is quite robust. 



3.1. CMB Dipole 

The dipole anisotropy stands apart from the rest of the CMB signal due to its large 
amplitude and to the understanding that it arises from our peculiar motion with respect to 
the CMB rest frame. In this section we present CMB dipole results based on a new analysis 
of the 5 year sky maps. Aside from an absolute calibration uncertainty of 0.2% (see §1]), the 
dominint source of uncertainty in the dipole estimate arises from uncertainties in Galactic 
foreground subtraction. Here we present results for two different removal methods: template- 
based cleaning an d an internal linear combination (ILC) of the WMAP multifrequency data 



( iGold et al.l 120081 ). Our final results are based on a combination of these methods with 



uncertainties that encompass both approaches. 

With template-based foreground removal, we can form cleaned maps for each of the 8 
high frequency DA's, Q1-W4, while the ILC method produces one cleaned map from a linear 
combination of all the WMAP frequency bands. We analyze the residual dipole moment in 
each of these maps (a nominal dipole based on the 3 year data is subtracted from the time- 
ordered data prior to map-making) using a Gibbs sampling technique which generates an 
ensemble of full-sky CMB realizations that are consistent with the data, as detailed below. 
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We evaluate the dipole moment of each full-sky realization and compute uncertainties from 
the scatter of the realizations. 

We prepared the data for the Gibbs analysis as follows. The N s ^ e = 512, template- 
cleaned maps were zeroed within the KQ85 mask, smoothed with a 10° FWHM Gaussian 
kernel, and degraded to iV^e = 16. Zeroing the masked region prior to smoothing prevents 
residual cleaning errors within the mask from contaminating the unmasked data. We add 
random white noise (12 /xK rms per pixel) to each map to regularize the pixel-pixel covariance 
matrix. The = 512 ILC map was also smoothed with a 10° FWHM Gaussian kernel 
and degraded to A^ si d e = 16, but the data within the sky mask were not zeroed prior to 
smoothing. We add white noise of 6 fiK per pixel to the smoothed ILC map to regularize 
its covariance matrix. Note that smoothing the data with a 10° kernel reduces the residual 
dipole in the maps by ~0.5%. We ignore this effect since the residual dipole is only ~0.3% 
of the full dipole amplitude to start with. 

The Gibbs sampler was run for 10,000 steps for each of the 8 template-cleaned maps 
(Q1-W4) and for each of 6 independent noise realizations added to the ILC map. In both 
cases we applied the KQ85 mask to the analysis and truncated the CMB power at / max = 32. 
The resulting ensembles of 80,000 and 60,000 dipole samples were analyzed independently 
and jointly. The results of this analysis are given in Table HI The first row combines the 
results from the template-cleaned DA maps; the scatter among the 8 DA's was well within 
the noise scatter for each DA, so the Gibbs samples for all 8 DA's were combined for this 
analysis. The results for the ILC map are shown in the second row. The two methods 
give reasonably consistent results, however, the Galactic longitude of the two dipole axis 
estimates differ from each other by about 2a. Since we cannot reliably identify one cleaning 
method to be superior to the other, we have merged the Gibbs samples from both methods to 
produce the conservative estimate shown in the bottom row. This approach, which enlarges 
the uncertainty to emcompass both estimates, gives 

(d, I, b) = (3.355 ± 0.008 mK, 263.99° ± 0.14°, 48.26° ± 0.03°), (1) 

where the amplitude estimate includes t he 0.2% absolute ca libration uncertainty. Given the 



CMB monopole temperature of 2.725 K (IMather et al.lll999l ). this amplitude implies a Solar 



System peculiar velocity of 369.0 ± 0.9 km s 1 with respect to the CMB rest frame. 



4. CALIBRATION IMPROVEMENTS 



With the 5 year processing we have refined our procedure for evaluating the instrument 
calibration, and have improved our estimates for the calibration uncertainty. The funda- 
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mental calibration source is sti ll the dipole anisotropy induced by WM AP's motion with 
respect to the CMB rest frame (IHinshaw et al.ll2003l ; iJarosik et al.ll2007l ). but several details 



of the calibration fitting have been modified. The new calibration solution is consistent with 
previous results in the overlapping time range. We estimate the uncertainty in the absolute 
calibration is now 0.2% per differencing assembly. 

The basic calibration procedure posits that a single channel of time-ordered data, di, 
may be modeled as 

dt = 9i [AT m + AT ai ] + b h (2) 

where % is a time index, and 6j are the instrument gain and baseline, at time step i, AT vi 
is the differential dipole anisotropy induced by WMAP's motion, and AT a i is the differential 
sky anisotropy. We assume that AT vi is known exactly and has the form 

T 



AT 



-Vj • [(1 + x im )n A , 



Xr 



(3) 



where T = 2.725 K is the CMB temperature (IMather et al.lll999l ). c is the speed of light, 
Vj is WMAP's velocit y with respect to th e CMB rest frame at time step i, x im is the loss 
imbalance parameter (IJarosik et al.l 120071 ). and n^j, and yib,% are the unit vectors of the 
A- and B-side lines of sight at time step i (in the same frame as the velocity vector). The 
velocity may be decomposed as 



v i — v WMAP-SSB,i + V SSB-CMB, 



(4) 



where the first term is WMAP's velocity with respect to the solar system barycenter, and the 
second is the barycenter velocity with respect to the CMB. The former is well determined 
from ep hemeris data, while the latter has been measured by COBE-DMR with an uncertainty 
of 0.7% (IKogut et al.lll996l ). Since the latter velocity is constant over WMAP's life span, any 
error in our assumed value of vssb-cmb will, in theory, be absorbed into a dipole contribution 
to the anisotropy map, T a . We test this hypothesis below. The differential sky signal has 
the form 



AT ai = (1 + x im )[I a (p Aji ) + P a (pA,i, lA,i 



X im )[I a (PB,i) + P a {PB,i, Tb,*)], (5) 



where p A ,i is the pixel observed by the A-side at time step i (and similarly for B), I a (p) is 
the temperature anisotropy in pixel p (the intensity Stokes parameter, 7), and P a (p, 7) is 
the polarization anisotropy in pixel p at polarization angle 7 (IHinshaw et al.ll2003l ) which is 
related to the linear Stokes parameters Q and U by 

P a (p, 7) = Q(p) cos 2 7 + U(p) sin 2 7 . (6) 
We further note that, in general, I a and P a depend on frequency owing to Galactic emission. 
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A main goal of the data processing is to simultaneously fit for the calibration and sky 
signal. Unfortunately, since the data model is nonlinear and the number of parameters 
is large, the general problem is intractable. In practice, we proceed iteratively as follows. 
Initially we assume the gain and baseline are constant for a given time interval, typically 
between 1 and 24 hours, 



9% — G k T~k < U < Tfc+i 

h = B k T k < ti < T k+X , 



(7) 



where £, is the time of the zth individual observation, and r k is the start time of the kth 
calibration interval. Throughout the fit we fix the velocity-induced signal, equation (J3J), 
using v SSB -cmb = [—26.29, —244.96, +275.93] km s _1 (in Galactic coordinates), and, for the 
first iteration, we assume no anisotropy signal, AT a = 0. Then, for each calibration interval 
k we perform a linear fit for G k and B k with fixed AT V + AT a . As we proceed through 
the intervals, we apply this calibration t o the raw data and ac cumulate a new estimate of 
the anisotropy map as per equation 19 of iHinshaw et al.l (120031 ). The procedure is repeated 
with each updated estimate of AT a . Once the calibration solution has converged, we fit the 
gain data, G k , to a model that is parameterized by the instrume nt detector vol t age a nd the 
temperatures of the receiver's warm and cold stages, equation 2 of lJarosik et al.l (120071 ). This 
parametrization still provides a good fit to the G k data, so we have not updated its form for 
the 5 year a nalysis. The updat ed best-fit parameters are given in the 5 year Explanatory 
Supplement (ILimon et al.ll2008l ). Note that for each radiometer, the relative gain vs. time 
over 5 years is determined by just two parameters. 

For the 5 year processing we have focused on the veracity of the "raw" calibration, G k 
and B k . Specifically, we have improved and/or critically reexamined several aspects of the 
iterative fitting procedure: 



We have incorporated the effect of far sidelobe picku p directly into the iterative cali- 
bration procedure, rather than as a fixed correction (jJarosik et al.l 120071 ). We do this 
by segregating the differential signal into a main beam contribution and a sidelobe 
contribution, 



ATi AT^^in i -\- AT'sifje j. 



(9) 



(IHill et al.ll2008l discuss how this segregation is defined in the 5 year processing.) After 
each iteration of the calibration and sky map estimation, we (re) compute a database 
of AT side on a grid of pointings using the new estimate of I a . We then interpolate the 
database to estimate AT side)i for each time step i. Note that AT side includes contri- 
butions from both the velocity-induced signal and the intrinsic anisotropy. Ignoring 
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sidelobe pickup can induce gain errors of up to 1.5% in K band, 0.4% in Ka band, and 
-0.25% in Q-W bands. 



In ge neral, the different channels within a DA have different center frequencies ( IJarosik et al 



20031 ) ; hence the different channels measure a slightly different anisotropy signal due 
to differences in the Galactic signal. We assess the importance of accounting for this 
in the calibration procedure. 

A single DA channel is only sensitive to a single linear polarization state. (WMAP mea- 
sures polarization by differencing orthogonal polarization channels.) Thus we cannot 
reliably solve for both P a and for I a at each channel's center frequency. We assess the 
relative importance of accounting for one or the other on both the gain and baseline 
solutions. 

We examine the sensitivity of the calibration solution to the choice of vssb-cmb and 
to assumptions of time-dependence in the gain. 



4.1. Calibration Tests 

We use a variety of end-to-end simulations to assess and control the systematic effects 
noted above. We summarize a number of the key tests in the remainder of this section. 

The first case we consider is a noiseless simulation in which we generate time-ordered 
data from an input anisotropy map which includes CMB and Galactic foreground signal (one 
map per channel, evaluated at the center frequency of each channel) and a known dipole 
amplitude. The input gain for each channel is fixed to be constant in time. We run the 
iterative calibration and sky map solver allowing for an independent sky map solution at 
each channel (but no polarization signal). When fitting for the calibration, we assume that 
vssb-cmb differs from the input value by 1% to see if the known, modulated velocity term, 
v wmap-ssb> properly "anchors" the absolute gain solution. The results are shown in the top 
panel of Figure [6] where it is shown that the absolute gain recovery is robust to errors in 
v ssb-cmb- We recover the input gain to better than 0.1% in this instance. 

The second case we consider is again a noiseless simulation that now includes only dipole 
signal (with Earth- v elocity modulation ), but here we vary the input gain using the flight- 



derived gain model (IJarosik et al.l 120071 ). The iterative solver was run on the K band data 
for 1400 iterations, again starting with an initial guess that was in error by 1%. The results 
are shown in the bottom panel of Figure [61 which indicate systematic convergence errors of 
>0.3% in the fitted amplitude of the recovered gain model. Since the input sky signal in this 
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case does not have any Galactic foreground or polarization components, we cannot ascribe 
the recovery errors to the improper handling of those effects in the iterative solver. We have 
also run numerous other simulations that included various combinations of instrument noise, 
CMB anisotropy, Galactic foreground signal (with or without individual center frequencies 
per channel), polarization signal, and input gain variations. The combination of runs are 
too numerous to report on in detail, and the results are not especially enlightening. The 
most pertinent trend we can identify is that when the input value of v SSB -cmb is assumed 
in the iterative solver, the recovered gain is in good agreement with the input, but when the 
initial guess is in error by 1%, the recovered gain will have comparable errors. We believe 
the lack of convergence is due to a weak degeneracy between gain variations and the sky map 
solution. Such a degeneracy is difficult to diagnose in the context of this iterative solver, 
especially given the computational demands of the system, so we are assessing the system 
more directly with a low-resolution parameterization of the gain and sky signal, as outlined 
in Appendix lAl 

Since the latter effort is still underway, we have adopted a more pragmatic approach 
to evaluating the absolute gain and its uncertainty for the 5 year data release. We proceed 
as follows: after 50 iterations of the calibration and sky map solver, the dominant errors 
in the gain and sky map solution are 1) a dipole in the sky map, and 2) a characteristic 
wave form that reflects a relative error between vssb-cmb and vwmap-ssb- At this point 
we can calibrate the amplitude of the gain error wave form to the magnitude of the velocity 
error in yssR-nivm- W e can then fit the gain solution to a linear combination of the gain 



model of iJarosik et al.l (120071 ) and the velocity error wave form. See Appendix |B] for details 
on this fitting procedure. In practice this fit is performed simultaneously on both channels 
of a radiometer since those channels share one gain model parameter. We have tested this 
procedure on a complete flight-like simulations that includes every important effect known, 
including input gain variations. The results of the gain recovery are shown in Figure [7J 
and based on this we conservatively assign an absolute calibration uncertainty of 0.2% per 
channel for the 5 year WMAP archive. 



4.2. Summary 

The series of steps taken to arrive at the final 5 year calibration are as follows: 

• Run the iterative calibration and sky map solver over the full 5 year data set for 50 
iterations, using 24 hour calibration intervals. This run starts with I a = P a = and 
updates I a for each individual channel of data. P a is assumed to be throughout this 
run. We keep the gain solution, G^, from this run and discard the baseline solution. 
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Run the iterative calibration and sky map solver over the full 5 year data set for 50 
iterations, using 1 hour calibration intervals. This run starts with I a = P a = and 
updates both using th e intensity and polariz ation data in the two radiometers per DA, 



as per Appendix D of iHinshaw et all (120031 ). We keep the baseline solution, B k , from 



this run and discard the gain solution. Both of these runs incorporate the sidelobe 
correction as noted above. 

Fit the gain solution, Gk simultaneously for the gain model and for an error in the ve- 
locity, Avssb-cmb, as described in Appendix [B] This fit is performed on two channels 
per radiometer with the gain model parameter T common to both channels. 

We average the best-fit velocity error over all channels within a frequency band under 
the assumption that the dipole is the same in each of these channels. We then fix the 
velocity error to a single value per frequency band and re-fit the gain model parameters 
for each pair of radiometer channels. 



Based on end-to-end simulations with flight-like noise, we estimate the absolute gain 
error per radiometer to be 0.2%. We believe the limiting factor in this estimate is a weak 
degeneracy between thermal variations in the instrument gain, which are annually modu- 
lated, and annual variations induced by errors in vssb-cmb- Since there is a small monotonic 
increase in the spacecraft temperature, additional years of data should allow improvements 
in our ability to separate these effects. 

Once we have finalized the gain model, we form a calibrated time-ordered data archive 
using the gain model and the 1 hour baseline estimates to calibrate the data. This archive 
also has a final estimate of the far sidelobe pickup subtracted from each time-ordered data 
point. However, we opt not to subtract a dipole estimate from the archive at this stage in 
the processing. 



5. BEAM IMPROVEMENTS 

In addition to reassessing the calibration, the other major effort undertaken to improve 
the 5 year data processing was to extend the physical optics model of the WMAP telescope 



based on flight measurements of Jupiter. This work is described in detail in lHill et al.l (120081 ) 
so we only summarize the key results with an emphasis on their scientific implications. The 
basic aim of the work is to use the flight beam maps from all 10 DA's to determine the 
in-flight distortion of the mirrors. This program was begun for the A-side mirror during the 
3 year analysis; for the 5 year analysis we have quadrupled the number of distortion modes 
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we fit (probing distortion scales that are half the previous size), and we have developed a 
completely new and independent model of the B-side distortions, rather than assuming that 
they mirror the A-side distortions. We have also placed limits on smaller scale distortions 
by comparing the predicted beam response at large angles to sidelobe data collected during 
WMAP's early observations of the Moon. 

Given the best-fit mirror model, we compute the model beam response for each DA 
and use it in conjunction with the flight data to constrain the faint tails of the beams, 
beyond ~ 1° from the beam peak. These tails are difficult to constrain with flight data 
alone because the Jupiter signal to noise ratio is low, but, due to their large areal extent 
they contain a non- negligible fraction (up to 1%) of the total beam solid angle. An accurate 
determination of the beam tail is required to properly measure the ratio of sub-degree-scale 
power to larger-scale power in the diffuse CMB emission (and to accurately assign point 
source flux). 



Figure 14 in iHill et all (120081 ) compares the beam radial profiles used in the 3 year and 
5 year analyses, while Figure 13 compares the /-space transfer functions derived from the 
Legendre transform of the radial profile. The important changes to note are the following. 

1. In both analyses we split the beam response into main beam and far sidelobe contri- 
bution s. In the 5 year analysis we have enlarged the radius at which this transition is 



made (iHill et al.ll2008l ). In both cases, we correct the time-ordered data for far sidelobe 
pickup prior to making sky maps, while the main beam contribution is only accounted 
for in the analysis of sky maps, e.g., in power spectrum deconvolution. As a result, 
the sky maps have a slightly different effective resolution which is most apparent in 
K band, as in Figure [51 However, in each analysis, the derived transfer functions are 
appropriate for the corresponding sky maps. 

2. In the 3 year analysis, the main beam profile was described by a Hermite polynomial 
expansion fit to the observations of Jupiter in the time-ordered data. This approach 
was numerically problematic in the 5 year analysis due to the larger transition radius; 
as a result, we now simply co-add the time-ordered data into radial bins to obtain the 
profiles. In both cases, the underlying time-ordered data is a hybrid archive consisting 
of flight data for points where the beam model p redicts a value a bove a given contour, 



and model values for points below the contour (iHill et al.ll2008l ). With the improved 
beam models and a new error analysis, we have adjusted these hybrid contours down 
slightly, with the result that we use proportionately more flight data (per year) in the 
new analysis. The radius at which the 5 year profile becomes mod el dominat e d (>5 0% 



of the points in a bin) is indicated by dotted lines in Figure 14 of IHill et ajj (120081 ). 
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The right column of Figure 14 in iHill et all (120081 ) shows the fractional change in solid 
angle due to the updated profiles. The main point to note is the ~1% increase in the 
V2 and W band channels, primaril y arising in the b in from 1 to 2 degrees off the beam 
peak. As can be seen in Figure 3 of IHill et all (120081 ). this is the angular range in which 
the new beam models produced the most change, owing to the incorporation of smaller 
distortion modes in the mirror model. The 3 year analysis made use of the model in 
this angular range which, in hindsight, was suppressing up to ~1% of the solid angle 
in the V and W band beams. (The longer wavelength channels are less sensitive to 
distortions in this range, so the change in solid angle is smaller for K-Q bands.) In the 
5 year analysis, we use relatively more flight data i n this regime, so we are less sensitive 
to any remaining model uncertainties. IHill et al.1 (120081 ) place limits on residual model 
errors and propagate those errors into the overall beam uncertainty. 



4. Figure 13 in IHill et al.1 (120081 ) compares the beam transfer functions, bi, derived by 
transforming the 3 year and 5 year radial profiles. (To factor out the effect of changing 
the transition radius, the 3 year profiles were extended to the 5 year radius using the 
far sidelobe data, for this comparison.) Since the transfer functions are normalized to 
1 at I — 1, the change is restricted to high I. In V and W bands, b\ has decreased by 
~0.5 - 1% due largely to the additional solid angle picked up in the 1-2 degree range. 
This amounts to a ~1 o change in the functions, as indicated by the red curves in the 
Figure. 



The calibrated angular power spectrum is proportional to l/g 2 bf, where g is the mean 
gain and b\ is the beam transfer function, thus the net effect of the change in gain and beam 
d etermination s is to increase the power spectrum by ~0.5% at I < 100, and by ~2.5% at high 



I. No 



ta et al. (2008) give a detai led evaluation of the power spectrum while iDunkley et al.l 



(120081 ) and iKomatsu et al.l (120081 ) discuss the implications for cosmology. 



LOW-/ POLARIZATION TESTS 



The 3 year data release included the first measurement of microwave polarization over 
the full sky, in the form of Stokes Q and U maps in each of 5 bands. The analysis of WMAP 
polarization data is complicated by the fact that the instrument was not designed to be a 
true polarimeter, thus a number of systematic effects had to be understood prior to assigning 
reliable error estimates to the data. iPage et al.1 (120071 ) presented the 3 year polarization data 
in great detail. In this section we extend that analysis by considering some additional tests 
that were not covered in the 3 year analysis. We note that all of the tests described in 
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this section have been performed on the template-cleaned reduced-foreground maps except 
for the final test of the Ka band data, described at the end of the section, which tests an 
alternative cleaning method. 



6.1. Year-to-Year Consistency Tests 

With 5 years of data it is now possible to subject the data to more stringent consistency 
tests than was previously possible. In general, the number of independent cross-power spectra 
we can form within a band with Nd differencing assemblies is Nd(Nd — l)/2 x N y + iVj x 
Ny(N y — l)/2. With 5 years of data, this gives 10 independent estimates each in K and Ka 
band, 45 each in Q and V band, and 190 in W band. For cross power spectra of distinct 
band pairs, with Ndi and Nd2 DA's in each band, the number is NdiNd2 x Ny. This gives 
50 each in KaQ and KaV, 100 each in KaW and QV, and 200 each in QW and VW. (For 
comparison, the corresponding numbers are 3, 15, & 66, and 18, 36, & 72 with 3 years of 
data.) 

We have evaluated these individual spectra from the 5 year data and have assign e d nois e 



uncertainties to each estimate using the Fisher formalism described in iPage et all (120071 ). 
We subject the ensemble to an internal consistency test by computing the reduced \ 2 of the 
data at each multipole I within each band or band pair, under the hypothesis that the data 
at each multipole and band measures the same number from DA to DA and year to year. 
The results of this test are given in Table [5] for the foreground-cleaned EE, EB, and BB 
spectra from I — 2 — 10 for all band pairs from KaKa to WW. There are several points to 
note in these results. 



1. For I > 6, the most significant deviation from 1 in reduced \ 2 i m an y spectrum or 
band, is 1.594 in the / = 7 BB spectrum for KaQ. With 50 degrees of freedom, this 
is a 3 er deviation, but given that we have 150 / > 6 samples in the table, we expect 
of order 1 such value. Thus we conclude that the Fisher-based errors provide a good 
description of the DA-to-DA and year-to-year scatter in the I > 6 polarization data. If 
anything, there is a slight tendency to overestimate the uncertainties at higher I. 

2. For I < 5, we find 37 out of 120 points where the reduced x 2 deviates from 1 at more 
than 4 a significance, indicating excessive internal scatter in the data relative to the 
Fisher errors. However, all but 5 of these occur in cross-power spectra in which one 
or both of the bands contain W band data. If we exclude combinations with W band, 
the remaining 72 points have a mode in the reduced \ 2 distribution of 1 with a slight 
positive skewness due to the 5 points noted above, which all contain Q band data. This 
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may be a sign of slight foreground residuals contributing additional noise to the Q band 
data, though we do not see similar evidence in the Ka band spectra which would be 
more foreground contaminated prior to cleaning. For Ka-V bands, we believe that the 
Fisher errors provide an adequate description of the scatter in this I < 5 polarization 
data, but we subject polarization sensitive cosmological parameter estimates, e.g., the 
optical depth, to additional scrutiny in §6.31 



3. Of special note is I = 3 BB which, as noted in iPage et al.l (120071 ) . is the power spec- 
trum mode that is least modulated in the WMAP time-ordered data. This mode is 
therefore quite sensitive to how the instrument baseline is estimated and removed and, 
in turn, to how the 1/f noise is modeled. In the accounting above, the / = 3 BB data 
have the highest internal scatter of any low-/ polarization mode. In particular, every 
combination that includes W band data is significantly discrepant; and the two most 
discrepant non-W band points are also estimates of I = 3 BB. We comment on the W 
band data further below, but note here that the final co-added BB spectrum (based 
on Ka, Q, and V band data) does not lead to a significant detection of tensor modes. 
However, we caution that any surprising scientific conclusions which rely heavily on 
the WMAP I = 3 BB data should be treated with caution. 



Based on the analysis presented above, we find the W band polarization data is still 
too unstable at low-/ to be reliably used for cosmological studies. We cite more specific 
phenomenology and consider some possible explanations in the remainder of this section. 

The 5 year co-added W band EE spectrum is shown in FiguresEl in the form of likelihood 
profiles from / = 2 — 7. At each multipole we show two curves: an estimate based on 
evaluati ng the likelihood multipole by multipole, and an estimate based on the pseudo-Q 
method (IPage et al.ll2007l ). The best-fit model EE spectrum, based on the combined Ka, Q, 
and V band data is indicated by the dashed lines in each panel. Both spectrum estimates 
show excess power relative to the model spectrum, with the most puzzling multipole being 
/ = 7 which, as shown in Table has an internal reduced x 2 °f 1-015, for 190 degrees of 
freedom. This data has the hallmark of a sky signal, but that hypothesis is implausible for 
a variety of reasons (IPage et al.l 120071 ). It is more likely due to a systematic effect that is 
common to a majority of the W band channels over a majority of the 5 years of data. We 
explore and rule out one previously neglected effect in §6.21 It is worth recalling that / = 7 
EE, like / = 3 BB , is a mode that is relatively poorly measured by WMAP, as discussed in 



Page et al.l (120071 ); see especially Figure 16 and its related discussion. 



The W band BB data also exhibit unusual behavior at / = 2, 3. In this case, these two 
multipoles have internal reduced x 2 greater than 6, and the co-added / = 2 point is nearly 
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10 a from zero. However, with 190 points in each 5 year co-added estimate it is now possible 
to look for trends within the data that were relatively obscure with only 3 years of data. In 
particular, we note that in the I = 2 estimate, there are 28 points that are individually more 
than 5 a from zero and that all of them contain Wl data in one or both of the DA pairs in 
the cross power spectrum. Similarly for I = 3, there are 14 points greater than 5 a and all of 
those points contain W4 data in one or both of the DA pairs. We have yet to pinpoint the 
significance of this result, but we plan to study the noise properties of these DA's beyond 
what has been reported to date, and to sharpen the phenomenology with additional years 
of data. 



6.2. Emissivity Tests 

In this section we consider time dependent emission from the WMAP optics 
didate for explaining the excess W band "signal" seen in the EE spectrum, mostly at I = 7. 
In the end, the effect proved not to be significant, but it provides a useful illustration of a 
common-mode effect that we believe is still present in the W band polarization data. 

From a number of lines of reasoning, we know that the microwave emissivity of the 
mirrors is a few percent in W band, and that it scales with frequ ency roughly like v} ' 5 across 



the WMAP frequency range, as expected for a classical metal (IBorn Wolj|l980l ). Hence 



this mechanism has the potent ial to explain a comm on- mode effect that is primarily seen in 



W band. Further, Figure 1 in iJarosik et al.l (120071 ) shows that the physical temperature of 



the primary mirrors are modulated at the spin period by ~200 fiK, with a dependence on 
solar azimuth angle that is highly repeatable from year to year. We believe this modulation 
is driven by solar radiation diffracting around the WMAP sun shield reaching the tops of 
the primary mirrors, which are only a few degrees within the geometric shadow of the sun 
shield. In contrast, the secondary mirrors and feed horns are in deep shadow and show no 
measurable variation at the spin period, so that any emission they produce only contributes 
to an overall radiometer offset, and will not be further considered here. 

As a rough estimate, the spin modulated emission from the primary mirrors could 
produce as much as ~ 0.02 x 200 = 4 /iK of radiometric response in W band, but the 
actual signal depends on the relative phase of the A and B-side mirror variations and the 
polarization state of the emission. In more detail, the differential signal, d(t), measured by 
a radiometer with lossy elements is 

d(t) = (1 - e A ) T A (t) - (1 - e B ) T B (t) + e A T p A (t) - e B T p B (t) (10) 

where e A = e A + e A + e A i s the combined loss in the A-side optics: (p)rimary plus (s)econdary 
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mirrors, plus the (f)eed horn, and likewise for the B-side. Ta,b is the sky temperature in the 
direction of the A or B-side line-of-sight; and T\ B is the physical temperature of the A or 
B-side primary mirror. 

The first two terms are the sky signal attenuated by the overall loss in the A and B 
side optics, respectiv ely. The effects of loss im balance, which arise when €a 7^ e#, have been 



studied extensively (IJarosik et al.l 120031 . 120071 ). We account for loss imbalance in the data 



processing and we marginalize over residual uncerta inties in the imbala nce coefficients when 



we form the pixel-pixel inverse covariance matrices (IJarosik et al.l 120071 ). Updated estimates 



of the loss imbalance coefficients based on fits to the 5 year data are reported in Table [6j 

In the remainder of this section we focus on the last two emissive terms in Equation [TUJ 
Recall that a WMAP differencing assembly consists of two radiometers, 1 and 2, that are 
sensitive to orthogonal linear polarization modes. The temperature and polarization signals 
are extracted by forming the sum and difference of the two radiometer outputs; thus, the 
emission terms we need to evaluate are 

m ± di{ t ) = (%^) n - (%^) n (n) 

where e A1 is the A-side primary mirror emissivity measured by radiometer 1, and so forth. 
The factor of 1 — e in the denominator applies a small correction for the mean loss, e = 
(ca + £b)/2, and arises from the process of calibrating the data to a known sky brightness 
temperature (§H). Note that we only pick up a polarized response if e x 7^ e 2 . 

We have simulated this signal in the time-ordered data using the measured primary 
mirror temperatures as template inputs. The emissivity coefficients were initially chosen to 
be consistent with the loss imbalance constraints. However, in order to produce a measurable 
polarization signal, we had to boost the emissivity differences to the point where they became 
unphysical, that is |ei — e 2 | > \e\ + e 2 |. Nonetheless, it was instructive to analyze this 
simulation by binning the resulting data (which also includes sky signal and noise) as a 
function of solar azimuth. The results are shown in the top panel of Figure [9] which shows 3 
years of co-added W band polarization data, the d\ — d 2 channel; the input emissive signal 
is shown in red for comparison. We are clearly able to detect such a signal with this manner 
of binning. We also computed the low-/ polarization spectra and found that, despite the 
large spin modulated input signal, the signal induced in the power spectrum was less than 2 
fiK 2 in /(/ + l)C i E£ '/27r, which is insufficient to explain the I = 7 feature in the W band EE 
spectrum. 

In parallel with the simulation analysis, we have binned the flight radiometer data by 
solar azimuth angle to search for spin modulated features in the polarization data. The 
results for W band are shown in the bottom panel of Figure [9] for the 5 year data. While 
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the x 2 per degree of freedom relative to zero is slightly high, there is no compelling evidence 
for a coherent spin modulated signal at the ~2 fiK level. In contrast, the simulation yielded 
spin modulated signals of 5-10 fiK and still failed to produce a significant effect in the EE 
spectrum. Hence we conclude that thermal emission from the WMAP optics cannot explain 
the excess W band EE signal. In any event, we continue to monitor the spin modulated data 
for the emergence of a coherent signal. 



6.3. Ka Band Tests 

The analysis presented in §6.11 shows that the Ka band polarization data is comparable 
to the Q and V band data in its internal consistency. That analysis was performed on data 



that had been f oreground cleaned using the template method discussed in iPage et al.l (120071 ) 



and updated in iGold et al.l (120081 ). In order to assess whether or not this cleaned Ka band 
data is suitable for use in cosmological parameter estimation we subject it to two further 
tests: 1) a null test in which Ka band data is compared to the combined Q and V band data, 
and 2) a parameter estimation based solely on Ka band data. 

For the null test, we form polarization maps by taking differences, |<Ska — §<Sqv> where 
S = Q,U are the polarization Stokes parameters, Sku are the maps formed from the Ka band 
data, and Sqv are the maps formed from the optimal combination of the Q and V band 
data. We evaluate the EE power spectrum from these null maps by evaluating the likelihood 
mode by mode while holding the other multipoles fixed at zero. The results are shown in 
Figure [TDl along with the best-fit model spectrum based on the final 5 year ACDM analysis. 
The spectrum is clearly consistent with zero, but to get a better sense of the power of this 
test, we have also used these null maps to estimate the optical depth parameter, r. The 
result of that analysis is shown as the dashed curve in Figure [TT1 where we find that the 
null likelihood peaks at r = and excludes the most-likely cosmological value with ~95% 
confidence. 

As a separate test, we evaluate the r likelihood using only the template-cleaned Ka 
band signal maps. The result of that test is shown as the blue curve in Figure [111 While the 
uncertainty in the Ka band estimate is considerably larger than the combined QV estimate 
(shown in red), the estimates are highly consistent. The result of combining Ka, Q, and V 
band data is shown in the black curve. 



Dunkley et al.l (120081 ) present a complementary method of foreground cleaning that 
makes use of Ka band data, in conjunction with K, Q, and V band data. Using a full 
6 parameter likelihood evalutaion, they compare the optical depth inferred from the two 
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cleaning methods wh ile using the full combined data sets in both cases: see Figure 9 of 



Dunkley et al.l (120081 ) for details. Based on these tests, we conclude that the Ka band data 



is sufficiently free of systematic errors and residual foreground signals that it is suitable for 
cosmological studies. The use of this band significantly enhances the overall polarization 
sensitivity of WMAP. 



7. SUMMARY OF 5- YEAR SCIENCE RESULTS 



De tailed present a tions o f the scientifi c resul ts from the 5 year data a re giv en bvlGold et al 



(p008h . lWright et al.l (120081 ) . iNolta et all (120081 ) . iDunklev et all (120081 ). and Komatsu et al.l 
( 120081 ). Sta rting with t he 5 y ear temperature and polarization maps, with their improved 
calibration, iGold et al.l (120081 ) give a new Markov Chain Monte Carlo-based analysis of fore- 
ground emission in the data. T heir results are broa dly consistent with previous analyses 
by the WMAP team and others (jEriksen et al.ll2007d ). while providing some new results on 
the microwave spectra of bright sources in the Galactic plane that aren't well fit by simple 
power-law foreground models. Figure [12] shows the 5 year CMB map based on the internal 
linear combination (LLC) method of foreground removal. 



Wright et al.l (120081 ) give a comprehensive analysis of the extragalactic sources in the 5 



year data, including a new analysis of variability made possible by the multi-year coverage. 
The 5 year WMAP source catalog now contains 390 objects and is reasonably complete 
to a flux of 1 Jy away fr om the Galactic plane. The new analysis of the WMAP beam 
response (IHill et al.l 120081 ) has led to more precise estimates of the point source flux scale 
for all 5 WMAP frequency bands. This information is incorporated in the new source 
catalog ( Wright et al. I2008T). and is also used to provide new brightness estimates of Mars, 
Jupiter, and Saturn (IHill et al.ll2008l ). We find significant (and expected) variability in Mars 
and Saturn over the course of 5 years a nd use that information to provide a preliminary 
recalibration of a Mars brightness model (IWrightll2007l ). and to fit a simple model of Saturn's 
brightness CIS db function of ring inclination. 



The temperature and polarization power spectra are presented in INolta et al.l (120081 ). 
The spectra are all consistent with the 3 year results with improvements in sensitivity com- 
mensurate with the additional integration time. Further improvements in our understanding 
of the absolute calibration and beam response have allowed us to place tighter uncertainties 
on the power spectra, over and above the reductions from additional data. These changes 
are all reflected in the new version of the WMAP likelihood code. The most notable im- 
provements arise in the third acoustic peak of the TT spectrum, and in all of the polarization 
spectra; for example, we now see unambiguous evidence for a 2nd dip in the high-Z TE spec- 
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trum, which further constrains deviations from the standard ACDM model. The 5 year TT 
and TE spectra are shown in Figure [131 We have also generated new maximum likelihood 
estimates of the low-/ polarization spectra: TE, TB, EE, EB, an d BB to complement our 
earlier estimates based on pseudo-Q methods (INofta et al.ll2008l ). The TB, EB, and BB 
spectra remain consistent with zero. 

The cosmologica l im plications of the 5 ye ar WMAP data are discussed in detail in 
Dunkley et al.l (120081 ) and iKomatsu et al.l (120081 ) . The now-standard cosmological model: a 
flat universe dominated by vacuum energy and dark matter, seeded by nearly scale- invariant, 
adiabatic, Gaussian random-phase fluctuations, continues to fit the 5 year data. WMAP has 
now determined the key parameters of this model to high precision; a summary of the 5 year 
parameter results is given in Table [71 The most notable improvements are the measurements 
of the dark matter density, Q c h 2 , and the amplitude of matter fluctu ations today, a^. Th e 
former is determined with 6% uncertainty using WMAP data only (IDunkley et al.ll2008l ). 
and with 3% uncerta inty when WMAP data is combined with BAO and SNe constraints 
( IKomatsu et al.ll2008l ). The latter is measured to 5% with WMAP data, and to 3% when 
combined with other data. The redshift of reionization is z re ion = 11.0 ± 1.4, if the universe 
were reionized instantaneously. The 2 a lower limit is 2 re i on >8.2, and instantaneous reion- 
ization at z rc ion = 6 is rejected at 3.5 a. The WMAP data continue s to favor models with a 
tilted primordial spectrum, n s = 0.963+q;q}5. IDunkley et al.l (120081 ) discuss how the ACDM 
model continues to fit a host of other astronomical data as well. 

Moving beyond t he standard ACDM m odel, when WMAP data is combined with BAO 
and SNe observations (IKomatsu et al.ll2008l ). we find no evidence for running in the spectral 
index of scalar fluctuations, dn s /dlnk = —0.028 ± 0.020 (68% CL). The new limit on the 
tensor-to-scalar ratio is r < 0.22 (95% CL), and we obtain tight, simultaneous limits on the 
(constant) dark energy equation of state and the spatial curvature of the universe: —0.14 < 
1 + w < 0.12 (95% CL) and -0.0179 < Q k < 0.0081 (95% CL). The angular power spectrum 
now exhibits the signature of the cosmic neutrino background: the number of relativistic 
degrees of freedom, expressed in units of the effective number of neutrino species, is found 
to be A/gff = 4.4 ±1.5 (68% CL), consistent with the standard value of 3.04. Models with 
N e ff = are disfavored at >99.5% confidence. A summary of the key cosmological parameter 
values is given in Table [7J, where we provide estimates using WMAP data alone and WMAP 
data combined with BAO and SNe observations. A complete tabulation of all parameter 
values for each model and dataset combination we studied is available on LAMBDA. 

The new data also place more stringent limits on deviati ons from Gaussianity , par- 
ity violations, and the amplitude of isocurvature fluctuations (IKomatsu et al.l 120081 ) . For 
example, new limits on physically motivated primordial non-Gaussianity parameters are 
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-9 < f l N C L al < HI (95% CL) and -151 < f^f < 253 (95% CL) for the local and equilateral 
models, respectively. 



8. CONCLUSIONS 

We have presented an overview of the 5 year WMAP data and have highlighted the 
improvements we have made to the data processing and analysis since the 3 year results were 
presented. The most substantive improvements to the processing include a new method for 
establishing the absolute gain calibration (with reduced uncertainty), and a more complete 
analysis of the WMAP beam response made possible by additional data and a higher fidelity 
physical optics model. Numerous other processing changes are outlined in £j2j 

The 5 year sky maps are consistent with the 3 year maps and have noise levels that are 
\/5 times less than the single year maps. The new maps are compared to the 3 year maps 
in §21 The main changes to the angular power spectrum are as follows: at low multipoles 
(I < 100) the spectrum is ~0.5% higher than the 3 year spectrum (in power units) due to the 
new absolute gain determination. At higher multipoles it is increased b y ~2.5%, due to the 



new beam response profiles, as explained in §5] and in iHill et al.l (120081 ). These changes are 
consistent with the 3 year uncertainties when one accounts for both the 0.5% gain uncertainty 
(in temperature units) and the 3 year beam uncertainties, which were incorporated into the 
likelihood code. 

We have applied a number of new tests to the polarization data to check internal con- 
sistency and to look for new systematic effects in the W band data (jEl). As a resu lt of these 



tests, and of new analyses of polarized foreground emission ( iDunkley et al.ll2008l ). we have 
concluded that Ka band data can be used along with Q and V band data for cosmological 
analyses. However, we still find a number of features in the W band polarization data that 
preclude its use, except in the Galactic plane where the signal to noise is relatively high. We 
continue to investigate the causes of this and have identified new clues to follow up on in 
future studies ( §6.ip . 

Scientific results gathered from the suite of 5 year papers are summarized in §7J The 
highlights include smaller uncertainties in the optical depth, r, due to a combination of 
additional years of data and to the inclusion of Ka band polarization data: instantaneous 
reionization at z Teion = 6 is now rejected at 3.5 o. New evidence favoring a non-zero neutrino 
abundance at the epoch of last scattering, made possible by improved measurements of the 
third acoustic peak; and new limits on the nongaussian parameter Jnl, based on additional 
data and the application of a new, more optimal bispectrum estimator. The 5 year data 
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continue to favor a tilted primordial fluctuation spectrum, in the range n s ~ 0.96, but a 
purely scale invariant spectrum cannot be ruled out at >3 a confidence. 

The WMAP observatory continues to operate at L2 as designed, and the addition of 
two years of flight data has allowed us to make significant advances in characterizing the 
instrument. Additional data beyond 5 years will give us a better understanding of the 
instrument, especially with regards to the W band polarization data since the number of 
jackknife combinations scales like the square of the number of years of operation. If W band 
data can be incorporated into the EE power spectrum estimate, it would become possible to 
constrain a second reionization parameter and thereby further probe this important epoch 
in cosmology. The WMAP data continues to uphold the standard ACDM model but more 
data may reveal new surprises. 



9. DATA PRODUCTS 

All of the WMAP data is released to the research community for further analysis through 



the Legacy Archive for Microwave Background Data Analysis (LAMBDA) at http:/ /lambda.gsfc. nasa.gov 
The products include the complete 5 year time-ordered data archive (both raw and cal- 
ibrated); the calibrated sky maps in a variety of processing stages (single year by DA, 
multi-year by band, high resolution and low resolution, smoothed, foreground-subtracted, 
and so forth); the angular power spectra and cosmological model likelihood code; a full table 
of model parameter values for a variety of model and data sets (including the best-fit model 
spectra and Markov chains); and a host of ancillary data to support further analysis. The 
WMAP Explanatory Supplem ent provides detail ed information about the WMAP in-flight 



operations and data products (ILimon et al.ll2008l ). 
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A. FISHER MATRIX ANALYSIS OF CALIBRATION AND SKY MAP 

FITS 

A.l. Least Squares Calibration and Sky Model Fitting 

Let % be a time index in the time ordered data. Let gi be parameters for the gain, ai m 
be parameters for the temperature anisotropy and b k be parameters for the baseline offset. 

The model of the time-ordered data (TOD) is 

rm = g t [AT m + AT at ] + b t , (Al) 

where % is a time index, AT vi is the differential dipole signal at time step i, including the 
CMB dipole, and AT ai is the differential anisotropy signal at time step i. The parameters 
of the model are the hourly gain and baseline values, and the sky map pixel temperatures 
(which goes into forming AT a . We fit for them by minimizing 

X 2 = £^?A (A2) 

i a i 

where q is the raw data, in counts, and U{ is the rms of the ith observation, in counts. 
The Fisher matrix requires taking the second derivative of y 2 with respect to all parameters 
being fit. In order to reduce the dimensionality of the problem to something manageable, 
we expand the calibration and sky signal in terms of a small number of parameters. We can 



This preprint was prepared with the AAS IATgX macros v5.2. 
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write 



9i 
bi 

AT ni 



j 

= b k B ki , 



y^Q/m \Xlm{fl A i) - Y lm (h B i)} , 
Im 



(A3) 
(A4) 
(A5) 



where G and B are function of time (defined below), a\ m are the harmonic coefficients of the 
map, and fiAi is the unit vector of the A-side feed at time step i, and likewise for B. 

A reasonable set of basis functions for the gain and baseline allow for an annual mod- 
ulation and a small number of higher harmonics. Note that this does not include power at 
the spin or precession period, which might be an important extension to consider. For now 
we consider the trial set 



Gji 



1 j = 

cos jOi j = l,...,j max 

sin(j - j max )6»i j = j max + 1, . . . , 2j B 



(A6) 



and 



B 



ki 



where = tan l {n y /h x ) 
components are evaluated in ecliptic coordinates. 



1 k = 

cos k9i k — 1, . . . , A; max , (A7) 

sin (A; A; max )^j k k max -\- 1, . . . , 2/c max 

Here n is the unit vector from WMAP to the Sun, and the 



A. 2. Evaluation of the Fisher Matrix 

We wish to evaluate the 2nd derivative 

1 d 2 x 2 

2 dpidpj (A8) 
where Pi and pj are the parameters we are trying to fit. The needed first derivatives are 

1 d X 2 \- (Cj - m^Gfj [AT m + AT m ] 

2 dgi' ~ ^ a\ ' 1 J 
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1 d X 2 



Then 



2 dai' m > 

i ay 



-E 



(Cj - mj)ft [yj' m '(nAi) - Ij'm'(nBi)] 



2 dgi'dgi" 



E 



(AH) 



(A12) 



1 <9 2 X 2 = v^yi 

2 dgi'dai'm' 4^ 



ft [*j'm'(n,u) - ij'm'(nBi)] Gj/j [AT„j + AT a , 



+ 0^(^-7^) (A13) 
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1 d 2 x 2 v 

2 d^'W ^ 
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Bk'jGj'j [AT vi + AT a , 
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E 
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(A14) 
(A15) 
(A16) 
(A17) 



From this we can form the inverse covariance matrix 



c- 1 



/ 1 9 2 X 2 1_9V_ 1 d 2 X 2 \ 

1 2 dgi'dgi" 2 dgi'da lllm „ 2 dgi' db k " ' 



1_9V 



i ay i ay 

2<9a i / m /%?" 2 da v m ida v , m n 2 da Vm ,db k '' 

1 d 2 x 2 1 9 2 X 2 1JV. 

\ 2db k 'dgi" 2db k 'da lllm „ 2 db k ' db k " J 



(A18) 



where the gain and baseline blocks are (2j max + 1) x (2j max + 1), and the sky map block is 

(^max "I - 1) X (/max "I - 1) • 

If we decompose C _1 using SVD the parameter covariance matrix can be inverted to 
have the form 

C = E ~ V V ® ^« (A19) 

0/7 



where the are the singular values, and the are the columns of the orthogonal matrix 
V. In this form, the uncertainty in the linear combination of parameters defined by V^) is 

l/Wi. 
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B. CALIBRATION MODEL FITTING WITH GAIN ERROR TEMPLATES 

B.l. Gain Error From Calibration Dipole Error 

Consider a simple model where the input sky consists of only a pure fixed (CMB) dipole, 
described by the vector d c , and a dipole modulated by the motion of WMAP with respect 
to the Sun, described by the time-dependent vector d v (t). The raw data produced by an 
experiment observing this signal is 

c{ti) = g(ti)[At c (ti) + At v (U)] (Bl) 

where c(tj) is the TOD signal in counts, g(ti) is the true gain of the instrument and At m (ti) 
is the differential signal produced by each dipole component (m = c, v) at time given the 
instrument pointing at that time. Note that we have suppressed the explicit baseline and 
noise terms here for simplicity. 

Now suppose we calibrate the instrument using an erroneous CMB dipole, d' c = rd c = 
(1 + Ar)d c , where r is a number of order one (and Ar < 1 so we can ignore terms of order 
Ar 2 ). The fit gain, gf(t), will then roughly have the form 

9 ' (t) = id< + d„ w r» (t Wd„wr (B2) 

where the vertical bars indicate vector magnitude. Now define d = d c + d v and expand to 
1st order in Ar to get 

d(t) • d c 



9 fit) = 9(t) 



1 - Ar- 



(B3) 



d(f)-d(f). 

Note that the term (d • d c )/(d • d) is dominated by a constant component of order d 2 c /(d 2 c + 
d 2 ) ~ 0.99, followed by an annually modulated term that is suppressed by a factor of order 
d v /d c . Thus an erroneous calibration dipole induces a specific error in the fit gain that can 
be identified and corrected for, assuming the time dependence of the true gain is orthogonal 
to this form. 



B.2. Gain Model Fitting 

In theory, the way to do this is as follows. We have a set of data in the form of the 
fit gains, gfj for each calibration sequence i, and we have a gain model, G(t;p n ), which is a 
function of time and a set of model parameters p n . Ideally we would like to fit the model to 
the true gain, g(t), but since we don't know the true gain, the next best thing is to modify 
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the gain model to have the same modulation form as the dipole gains have and to fit for this 
modulation simultaneously with the other gain model parameters. Thus x 2 takes the form 

2 \- \9i-Gi(p n )f \- [g f ,i - Gi(p n )(l - Ar/j.j)] 2 

x = 2^ ~2 = 2^ ~2 ' ( B4 ) 

i * i * 

where /<y = (d ■ d c ) /(d-d) evaluated at time ij, or is a function generated from simulations. 

Since the system is nonlinear, it must be minimized using a suitable nonlinear least 
squares routine. However, we can analyze the parameter covariance matrix directly by 
explicitly evaluating the 2nd derivative of x 2 with respect to the model parameters 

1 <9V 

< B5 > 

First compile the necessary 1st derivatives 

1 dx 2 [9f,i - Gi(p n )(l - Arf dji )} (Gi fo) 

2dAr = ^ V 2 (B6) 

i 1 

1 d X 2 \- \m - Gi(Pn)(l - Arf dii )] (-dGi/dp m )(l - Arf d>i ) 

~2dp- = ^ V 2 (B7) 

(We evaluate the individual dG/dp m terms below.) Next the various 2nd derivatives are 

1 d 2 X 2 (Gif^idU 



2 dArdAr a 



Vfj jd t i)y^i jd,i) ^gg-j 



1 d 2 X 2 ^ (GiUM-dGj/dpmXl-ArUi) ^ _ 
9 7)A^T = ^ + ° U* ~ (B9) 



2 dArdp m *— J of 

1 d\ 2 ^ (dGi/dpmHl-ArUMdGj/dpnHl-Arfoi) ^ 
2dp~di n = \ of +0^- Gi ). (BIO) 

In the last two expressions, we neglect the term proportional to d 2 G / dp m dp n because the 
prefactor of — Gi) is statistically zero for the least squares solution. 



Finally, we evaluate the dG/dp m terms. The gain model has the form (jjarosik et al. 
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Gi = a — — — Q , (Bll) 

J-FPAKH) — J FPA 

where Tr XB = 290 K, and a, Vq, and Tp PA are parameters to be fit. The necessary 1st 
derivatives are 

dGt/da = ^'^f"' 1 "^ , (B12) 

J-FPA\ti) ~ J FPA 
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Fig. 1. — Five-year temperature sky maps in Galactic coordinates smoothed with a 0.2° 
Gaussian beam, shown in Mollweide projection, top: K band (23 GHz), middle-left: Ka 
band (33 GHz), bottom-left: Q band (41 GHz), middle-right: V band (61 GHz), bottom- 
right: W band (94 GHz). 
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Fig. 2. — Five-year Stokes Q polarization sky maps in Galactic coordinates smoothed to 
an effective Gaussian beam of 2.0°, shown in Mollweide projection, top: K band (23 GHz), 
middle-left: Ka band (33 GHz), bottom-left: Q band (41 GHz), middle-right: V band (61 
GHz), bottom-right: W band (94 GHz). 
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Fig. 3. — Five-year Stokes U polarization sky maps in Galactic coordinates smoothed to 
an effective Gaussian beam of 2.0°, shown in Mollweide projection, top: K band (23 GHz), 
middle-left: Ka band (33 GHz), bottom-left: Q band (41 GHz), middle-right: V band (61 
GHz), bottom-right: W band (94 GHz). 
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Fig. 4. — Five-year polarization sky maps in Galactic coordinates smoothed to an effective 
Gaussian beam of 2.0°, shown in Mollweide projection. The color scale indicates polarized 
intensity, P = \J Q 2 + U 2 , and the line segments indicate polarization direction in pixels 
whose signal-to-noise exceeds 1. top: K band (23 GHz), middle-left: Ka band (33 GHz), 
bottom-left: Q band (41 GHz), middle-right: V band (61 GHz), bottom-right: W band (94 
GHz). 
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Fig. 5. — Difference between the 5 year and 3 year temperature maps, left column: the 
difference in the maps, as delivered, save for the subtraction of a relative offset (Table [3]), 
right column: the difference after correcting the 3 year maps by a scale factor that accounts 
for the mean gain change, ~ 0.3%, between the 3 year and 5 year estimates, top to bottom: 
K, Ka, Q, V, W band. The differences before recalibration are dominated by galactic plane 
emission and a dipole residual: see Table [3j which also gives the changes for I — 2, 3. 
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Fig. 6. — Gain convergence tests using the iterative sky map & calibration solver run on a 
pair of simulations with known, but different, inputs. Both panels show the recovered gain 
as a function of iteration number for a 4-channel K band simulation. The initial calibration 
guess was chosen to be in error by 1% to test convergence; the output solutions, extrapolated 
with an exponential fit, are printed in each panel, top: Results for a noiseless simulation 
that includes a dipole signal (with Earth-velocity modulation) plus CMB and foreground 
anisotropy (the former is evaluated at the center frequency of each channel) . The input gain 
was set to be constant in time. The extrapolated solutions agree with the input values to 
much better than 0.1%. bottom: Results for a noiseless simulation that includes only dipole 
signal (with Earth-velocity modulation) but no CMB or foreground signal. In this case the 
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Fig. 7. — Gain error recovery test from a flight-like simulation that includes every effect 
known to be important. Using the daily dipole gains recovered from the iterative sky map 
& calibration solver as input, the gain convergence error, shown here, is fit simultaneously 
with the gain model parameters, not shown, following the procedure outlined in Appendix iBl 
The red trace indicates the true gain error for each WMAP channel, based on the known 
input gain and the gain solution achieved by the iterative solver on its final iteration. The 
black trace shows the gain error recovered by the fit, averaged by frequency band. The 
channel-to-channel scatter within a band is <0.1%, though the mean of Ka band error is of 
order 0.1%. 
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Fig. 8. — W band EE power spectrum likelihood from I = 2 — 7 using two separate estimation 
methods: black: maximum likelihood and red: pseudo-C^. The vertical dashed lines indicate 
the best-fit model power spectrum based on fitting the combined Ka, Q, and V band data. 
The two spectrum estimates are consistent with each other, except at I — 3. The maximum 
likelihood estimates are wider because they include cosmic variance whereas the pseudo-C; 
estimates account for noise only. Both estimates show excess power in the W band data 
re lative to the best-fi t model, and to the combined KaQV band spectrum, shown in Figure 6 
of iNofta et al.l (120081 ). The extreme excess in the I = 7 pseudo-C; estimate is not so severe 
in the maximum likelihood, but both methods are still inconsistent with the best-fit model. 
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Fig. 9. — top: Simulated W band data with a large polarized thermal emission signal injected, 
binned by solar azimuth angle. The red trace shows the input waveform based on the flight 
mirror temperature profile and a model of the polarized emissivity. The black profile is the 
binned co-added data which follows the input signal very well. The thickness of the points 
represents the 1 a uncertainty due to white noise, bottom: Same as the top panel but for 
the 5 year flight data. The reduced \ 2 °f t ne binned data with respect to zero is 2.1 for 36 
degrees of freedom, but this does not account for 1/f noise, so the significance of this result 
requires further investigation. However, the much larger signal in the simulation did not 
produce an EE spectrum with features present in the flight W band EE spectrum, so the 
feature in the binned flight data cannot account for the excess I = 7 emission. 
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Fig. 10. — The EE power spectrum computed from the null sky maps, |«Sk a — f 5qV) where 
S = Q,U are the polarization Stokes parameters, and Sqv is the optimal com bination of the 
Q and V band data. The pink curve is the best-fit theoretical spectrum from iDunkley et al. 
( 120081 ) . The spectrum derived from the null maps is consistent with zero. 
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Fig. 11. — Estimates of the optical depth from a variety of data combinations. The dashed 
curve labeled Null uses the same null sky maps used in Figure [101 The optical depth 
obtained from Ka band data alone (blue) is consistent with independent estimates from the 
combined Q and V band data (red). The final 5 year analysis uses Ka, Q, and V band 
data combined (black). These estimates all use a 1-parameter likelihood estimation, holding 
other parameters fixed except for t he fluctuation amp litude, which is adjusted to fit the first 
acoustic peak in the TT spectrum (IPage et al.l 120071). The degen eracy between r and other 
ACDM parameters is small: see Figure 7 of iDunkley et al.l (120081 ). 
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Fig. 12. — The foreground- reduced Internal Linear Combination (ILC) map based on the 5 
year WMAP data. 
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Fig. 13. — The temperature (TT) and temperature-polarization correlation (TE) power 
spectra based on the 5 year WMAP data. The addition of 2 years of data provide more 
sensitive measurements of the third peak in TT and the high-/ TE spectrum, especially the 
second trough. 



Table 1. Differencing Assembly (DA) Properties 



DA 


A a 


z/ a 




#FWHM° 


^o(I) d 


^o(Q,U) d 










(mm) 


(GHz) 




(°) 


(mK) 


(mK) 


(GHz) 


(GHz) 


(GHz) 


Kl 


13.17 


22.77 


1.0135 


0.807 


1.436 


1.453 


22.47 


22.52 


22.78 


Kaf 


9.079 


33.02 


1.0285 


0.624 


1.470 


1.488 


32.71 


32.76 


33.02 


Ql 


7.342 


40.83 


1.0440 


0.480 


2.254 


2.278 


40.47 


40.53 


40.85 


Q2 


7.382 


40.61 


1.0435 


0.475 


2.141 


2.163 


40.27 


40.32 


40.62 


VI 


4.974 


60.27 


1.0980 


0.324 


3.314 


3.341 


59.65 


59.74 


60.29 


V2 


4.895 


61.24 


1.1010 


0.328 


2.953 


2.975 


60.60 


60.70 


61.27 


Wl 


3.207 


93.49 


1.2480 


0.213 


5.899 


5.929 


92.68 


92.82 


93.59 


W2 


3.191 


93.96 


1.2505 


0.196 


6.565 


6.602 


93.34 


93.44 


94.03 


W3 


3.226 


92.92 


1.2445 


0.196 


6.926 


6.964 


92.34 


92.44 


92.98 


W4 


3.197 


93.76 


1.2495 


0.210 


6.761 


6.800 


93.04 


93.17 


93.84 



a Effective wavelength and frequency for a thermodynamic spectrum. 

b Conversion from antenna temperature to thermodynamic temperature, AT = g(u)AT^. 

c Full-width-at-half-maximum from radial profile of A- and B-side average beams. Note: 
beams are not Gaussian. 

d Noise per observation for resolution 9 and 10 /, Q, & U maps, to ~0.1% uncertainty. 
a(p)=a N; h 1 J 2 (p). 

c Effective frequency for synchrotron (s), free-free (ff), and dust (d) emission, assuming 
spectral indices of j3 — —2.9, —2.1, +2.0, respectively, in antenna temperature units. 
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Table 2. Lost and Rejected Data 



Category 




K-band 


Ka-banc 


I Q-band 


V-band 


W-band 


Lost or incomplete telemetry(%) 




0.12 


0.12 


0.12 


0.12 


0.12 


Spacecraft anomalies(%) 




0.44 


0.46 


0.52 


0.44 


0.48 


Planned stationkeeping maneuvers(%) 


0.39 


0.39 


0.39 


0.39 


0.39 


Planet in beam (%) 




0.11 


0.11 


0.11 


0.11 


0.11 


Total lost or rejected (%) 




1.06 


1.08 


1.14 


1.06 


1.10 


Table 3. Change in 


low-/ Power from 3 


year Data 






Band 


I = o a 


i = r 


I = 2 b I 


= 3 b 








(fK) 


(JJ.K) 


(/zK 2 ) GuK 2 ) 






K 


9.3 


5.1 


4.1 


0.7 






Ka 


18.9 


2.1 


2.8 


0.2 






Q 


18.3 


0.4 


2.5 


0.5 






V 


14.4 


7.3 


1.2 


0.0 






w 


16.4 


3.5 


1.0 


0.0 







a / = 0, 1 - Amplitude in the difference 
map, outside the processing cut, in /iK. 

b l = 2, 3 - Power in the difference map, 
outside the processing cut, 1(1 + 1) Ci/2n, 
in /iK 2 . 
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Table 4. WMAP 5 year CMB Dipole Anisotropy a 



Cleaning 


4 b 


dy 


d z 


d c 


I 


b 


method 


(mK) 


(mK) 


(mK) 


(mK) 


n 


(°) 


Templates 


-0.229 ±0.003 


-2.225 ±0.003 


2.506 ±0.003 


3.359 ±0.008 


264.11 ±0.08 


48.25 ±0.1 


ILC 


-0.238 ±0.003 


-2.218 ±0.002 


2.501 ±0.001 


3.352 ±0.007 


263.87 ±0.07 


48.26 ±0.1 


Combined 


-0.233 ±0.005 


-2.222 ±0.004 


2.504 ±0.003 


3.355 ±0.008 


263.99 ±0.14 


48.26 ±0.1 



a The CMB dipole components for two different galactic cleaning methods are given in the first two ro 
The Gibbs samples from each set are combined in the last row to produce an estimate with conserval 
uncertainties that encompasses both cases. 

b The cartesian dipole components are given in Galactic coordinates. The quoted uncertainties reflect 
effects of noise and sky cut, for illustration. An absolute calibration uncertainty of 0.2% should be addec 
quadrature. 

c The spherical components of the dipole are given in Galactic coordinates. In this case the quoted un< 
tainty in the magnitude, d, includes the absolute calibration uncertainty. 
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Table 5. Polarization x 2 Consistency Tests a 



Multipole KaKa KaQ KaV KaW QQ QV QW VV VW WW 
(10) b (50) (50) (100) (45) (100) (200) (45) (200) (190) 



EE 




2 


0.727 


3 


1.373 


4 


1.561 


5 


0.914 


6 


1.003 


7 


0.600 


8 


1.578 


q 


n 7fio 

U. 1 uw 


i n 


fl 4Q4 


EB 




2 


0.900 


3 


0.719 


4 


0.746 


5 


1.161 


6 


0.475 


7 


1.014 


8 


0.849 


9 


0.743 


10 


0.413 


BB 




2 


2.038 


3 


0.756 


4 


1.058 


5 


1.221 


6 


0.379 


7 


1.925 


8 


0.804 


9 


0.320 



1.059 


1.019 


1.301 


0.994 


1.683 


1.355 


1.816 


1.341 


2.033 


1.313 


1.062 


1.275 


0.847 


0.688 


1.124 


0.671 


0.689 


0.936 


1.262 


1.337 


1.212 


n 7i n 


n 8Qi 


D 890 




U. c/c/U 


n Qi4 


1.297 


1.179 


2.074 


1.599 


0.651 


2.182 


1.702 


1.378 


1.777 


0.948 


0.945 


1.003 


1.183 


0.651 


0.687 


1.007 


0.829 


0.700 


0.897 


1.279 


0.861 


0.734 


1.007 


1.112 


1.003 


1.316 


0.859 


1.570 


1.244 


2.497 


0.868 


0.808 


1.817 


1.455 


1.522 


2.144 


1.659 


1.742 


2.036 


0.805 


0.483 


0.812 


1.594 


0.967 


1.332 


1.005 


0.999 


0.912 


0.489 


0.502 


0.450 



1.586 


0.690 


1.179 


1.092 


1.614 


1.325 


0.993 


1.126 


1.581 


1.631 


1.052 


1.155 


0.740 


0.856 


1.049 


0.936 


0.780 


0.864 


1.080 


0.763 


0.608 


fl 589 


n 79fi 


n 65i 

\J.\JO ± 






U.OUU 


1.006 


0.915 


2.126 


1.295 


0.986 


2.739 


1.926 


1.110 


1.435 


1.149 


1.232 


1.468 


0.829 


1.023 


0.814 


0.817 


0.759 


1.112 


0.681 


0.689 


0.955 


0.820 


0.798 


0.686 


0.722 


0.900 


0.693 


1.340 


1.219 


2.529 


3.027 


1.717 


3.496 


1.007 


0.905 


1.786 


0.889 


1.057 


1.271 


1.009 


0.861 


1.238 


1.074 


0.817 


0.928 


1.069 


0.782 


0.831 


0.884 


0.491 


0.729 



0.894 


1.078 


1.152 


1.005 


1.386 


1.519 


1.195 


1.596 


1.724 


0.589 


0.881 


1.252 


1.384 


1.168 


1.142 


0.900 


1.064 


1.015 


1.025 


0.871 


0.749 


n 7Qi 


n 89i 


fl 7Q5 


U.Ul u 




n Q/n 

U. cf^tO 


1.242 


2.085 


2.309 


1.095 


3.276 


3.157 


1.028 


1.279 


1.861 


0.699 


1.122 


1.516 


1.201 


1.136 


0.960 


0.616 


0.802 


1.233 


1.021 


0.954 


0.996 


0.882 


0.808 


0.824 


1.124 


0.836 


0.852 


0.694 


1.631 


9.195 


0.601 


2.545 


5.997 


0.752 


1.403 


1.984 


1.078 


1.660 


1.255 


0.800 


0.767 


0.955 


0.772 


0.994 


1.024 


0.997 


0.879 


0.943 


0.748 


0.664 


0.959 



Table 5 — Continued 



Multipole 


KaKa 


KaQ 


KaV 


KaW 


QQ 


QV 


QW 


vv 


vw 


WW 




(10) b 


(50) 


(50) 


(100) 


(45) 


(100) 


(200) 


(45) 


(200) 


(190) 


10 


1.181 


1.162 


1.028 


0.980 


1.218 


1.165 


0.951 


1.079 


0.621 


0.791 



a Table gives y 2 per degree of freedom of the independent spectrum estimates per multipole 
per band or band-pair, estimated from the template-cleaned maps. See text for details. 

b Second header row indicates the number of degrees of freedom in the reduced \ 2 f° r that 
spectrum. See text for details. 



Table 6. Loss Imbalance Coefficients 3, 



DA 


•£im,l 


Xim,2 




(%) 


(%) 


Kl 


0.012 


0.589 


Kal 


0.359 


0.148 


Ql 


-0.031 


0.412 


Q2 


0.691 


1.048 


VI 


0.041 


0.226 


V2 


0.404 


0.409 


Wl 


0.939 


0.128 


W2 


0.601 


1.140 


W3 


-0.009 


0.497 


W4 


2.615 


1.946 



a Loss imbalance is 
defined as x- 1Ta = (e^ — 
e B )/{e A + e B ). See 
§6.21 and iJarosik et al. 



(120071 ) for details. 
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Table 7. Cosmological Parameter Summary 



Description Symbol WMAP-only WMAP+BAO+SN 



Parameters for Standard ACDM Model a 



Age of universe 




13.69 ±0.13 Gyr 


13.72 ±0.12 Gyr 


Hubble constant 


U 


71. 7 km/s/Mpc 


70.5 ± 1.3 km/s/Mj 


Baryon density 




0.0441 ± 0.0030 


0.0456 ±0.0015 


Physical baryon density 


n b h 2 


0.02273 ±0.00062 


0.02267+^59 


Dark matter density 


n c 


0.214 ±0.027 


0.228 ±0.013 


Physical dark matter density 


n c h 2 


0.1099 ±0.0062 


0.1131 ±0.0034 


Dark energy density 


n A 


0.742 ±0.030 


0.726 ±0.015 


Curvature fluctuation amplitude, fco = 0.002 Mpc -1 b 


A 2 

TZ 


(2.41 ± 0.11) x 10~ 9 


(2.445 ± 0.096) x 10 


Fluctuation amplitude at 8/i _1 Mpc 


°8 


0.796 ±0.036 


0.812 ±0.026 


1(1 + l)C^/2n 


C220 


5756 ± 42 ^K 2 


5751±43 (1K 2 


Scalar spectral index 


n s 


963+ 014 
u.yoo_ 015 


0.960 ±0.013 


Rcdshift of matter-radiation equality 




31 76 +151 
01 ' °-150 


3253+J 9 . 


Angular diameter distance to matter-radiation eq. c 


dA(Zcq) 


14279±i|| M P C 


14200+^ M P C 


Rcdshift of decoupling 




1090.51 ±0.95 


1090.88 ±0.72 


Age at decoupling 




38008ltg£? yr 


376971+J^ yr 


Angular diameter distance to decoupling c - d 


d A {z«) 


14115±i|f Mpc 


14034+.^ Mpc 


Sound horizon at decoupling d 


r s {z*) 


146.8 ± 1.8 Mpc 


145.9+i;2 M P C 


Acoustic scale at decoupling d 


Ia(z#) 


302.08t°;g 


302.13 ±0.84 


Rcionization optical depth 


T 


0.087 ±0.017 


0.084 ±0.016 


Rcdshift of rcionization 


^reion 


11.0 ± 1.4 


10.9 ± 1.4 


Age at reionization 


^reion 




432+;$ Myr 


Parameters for Extended Models e 




Total density f 


^tot 


i-09»±g:Sgg 


1 00 50+ 0060 
i.uuou_ 0061 


Equation of state s 


W 


-i.06lg:S 


-0 9 92+ 061 


Tensor to scalar ratio, fc = 0.002 Mpc" 1 h < h 


r 


< 0.43 (95% CL) 


< 0.22 (95% CL) 


Running of spectral index, fc = 0.002 Mpc -1 h ' 1 


dn s /dh\k 


-0.037 ±0.028 


-0.028 ±0.020 



Table 7 — Continued 



Description 


Symbol 


WMAP-only 


WMAP+BAO+SN 


Neutrino density J 




< 0.014 (95% CL) 


< 0.0071 (95% CL) 


Neutrino mass J 




< 1.3 eV (95% CL) 


< 0.67 eV (95% CL) 


Number of light neutrino families k 


Neff 


> 2.3 (95% CL) 


4.4 ± 1.5 



a The parameters reported in the first sect i on as sume the 6 parameter ACDM model, 
first using WMAP d ata only |Dunklev et all 120081 ). then using WMAP+BAO+SN data 
(jKomatsu et alj|2008t ). 

b fc = 0.002 Mpc- 1 < — ► l eS w 30. 

c Comoving angular diameter distance. 

d l A {z*) = ttcIa(z*) r s (z*)~ l . 

e The parameters reported in the second section place limits on deviations from the ACDM 
model, first using WM AP data only jPunklev et al.ll2008h . then using WMAP+BAO+SN data 
( Komatsu et al. 2008h . A complete listing of all parameter values and uncertainties for each 
of the extended models studied is available on LAMBDA. 

f Allows non-zero curvature, fifc 7^ 0. 

s Allows mi/-1, but assumes w is constant. 

h Allows tensors modes but no running in scalar spectral index. 

'Allows running in scalar spectral index but no tensor modes. 

J Allows a massive neutrino component, £l u + 1 0. 

k Allows iV e ff number of relativistic species. The last column adds the HST prior to the other 
data sets. 



